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1. Introduction 

As is well known, at nonzero baryon chemical potential the complexity of the fermion deter- 
minant prohibits the use of importance sampling in lattice QCD. This makes the determination of 
the QCD phase diagram an outstanding open problem (a summary of what is possible despite this 
obstruction can be found in Ref. [jl]]). Stochastic quantization [|2|] does not rely on importance sam- 
pling: the configurations that dominate in the partition function are found by integrating complex 
Langevin equations in an enlarged phase space, after a complexification of the degrees of freedom 
[Jl, ||, ||]. During last year's Conference I reported [§] on work done in collaboration with Ion- 
Olimpiu Stamatescu on QCD with static quarks and related one-link models 0. In the one-link 
models we found excellent results for all values of the chemical potential, from zero all the way to 
saturation. Moreover, the sign problem was present but clearly not dangerous. The results in the 
heavy dense limit of QCD were preliminary but encouraging. 

These results triggered a number of intriguing questions. For instance, can complex Langevin 
dynamics handle 

• a severe sign problem? 

• the thermodynamic limit? 

• phase transitions? 

• the Silver Blaze problem ||? 

I will answer these questions positively in a theory that is much simpler than QCD: the rel- 
ativistic Bose gas at nonzero chemical potential ijjj. Even though it lacks the intricacies of a 
nonabelian gauge theory, it is a four-dimensional lattice field theory with a sign problem, a Sil- 
ver Blaze problem, and a phase transition. As an additional benefit, one can study the complex 



Langevin equations perturbatively, allowing detailed analytical insight into the algorithm [10]. Of 



course, there are many more questions to ask and work to address those is currently in progress. 

2. Relativistic Bose gas 

The lattice action is 



v=l 



(2.1) 



where d = 4 and we take m 2 > 0. It satisfies 5*(ju) = 5(— jU*), just as in QCD. In continuum no- 
tation, the nonderivative part of the Lagrangian density reads (m 2 — /x 2 ) |0| 2 + A|0| 4 : at tree level, 
a phase transition at pL c = m separates the vacuum (symmetric) and the high-density (symmetry- 
broken) phases. The /I -independence of bulk thermodynamic quantities as long as pL < ;U C , even 
though the action, the Boltzmann weight, and hence the field configurations depend explicitly on 
jU, is the Silver Blaze problem [^]. Equating the imaginary part of the action with zero yields the 
phase quenched theory: this is a theory with a /I -dependent mass parameter and no sign problem. 
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After decomposing = (0i +i<f^) /y/2 and complexifying fl — > + 0* (a = 1,2), it is 
straightforward to write down the complex Langevin equations for the (now) four real fields, 



djl 

de 



-Re 



8S 

S(f>a 



de 



-Im 



_55 



SHtf 



(2.2) 



where is the Langevin time and 7] fl is Gaussian noise with {f] a f]b) = 25^. It is also straightfor- 
ward to solve these equations numerically [0]. I have used m = X = 1 on lattices of size Q. = N , 
with A 7 " = 4,6,8, 10, and Langevin time step e = 5 x 10~ 5 and not encountered instabilities in the 
numerical integration. The phase quenched theory is solved using real Langevin dynamics (with 
^=0). 

The resulting density (n) = Qr l d\nZ/d}X is shown in Fig. |], in the full theory (left) and the 
phase quenched theory (right). The Silver Blaze feature in the full theory is immediately visible. In 
the thermodynamic limit the density is strictly zero as long as /I < /i c « 1.15. The phase quenched 
theory has no Silver Blaze region: instead the density increases as soon as pL is nonzero. The 
/i -independence in the full theory is a direct result of the complexity of the action and the sign 
problem. 
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Figure 1: Density in the full (left) and the phase quenched (right) theory as a function of chemical potential 
on lattices with size N A , with N = 4,6,8, 10 (m = A = 1). 



In this model the complex Langevin equations can be studied perturbatively [jlOJ. Ignoring 
interactions, the solution of these equations reads, in momentum space and in the Silver Blaze 
region only, 



C( ) = e ~ A, ° [^{B p e)^l p {0) + i^m{B p d)e ab ^ Ki 

+ [ dse- A "^cos[B p (e- S )}T] a>p (s), (2.3) 
Jo 

<t>l P (0) = e- A » d [cos(B p d)^ p (0)-ism(B p G)e ab ^ p (0)] 

-ij dse-^-'hmiBpid-s^Sab^pis). (2.4) 
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Here <j> a ' p (0) are the initial conditions, z a b is the totally antisymmetric tensor with £\\ = 1, and the 
coefficients A p , B p are 

3 

A p = m 2 + 4^sin y + 2(1 -cosh /icos/? 4 ) m 2 - fx 2 + p 2 + pj, (2.5) 

i=i ^ 

Bp = 2 shlh /X sin /?4 —> 2/J.P4. (2.6) 

The expressions after the arrows correspond to the formal continuum limit. This solution is only 
valid when A p > (or /I < pL L ), which reflects the standard instability of the free Bose gas. When 
/I > /i c , symmetry breaks and the inclusion of the A|0| 4 term is necessary for stabilization. 

From the analytical solution, we find that both the independence of initial conditions and 
convergence at large Langevin time are controlled by e~ Ap6 . This is demonstrated numerically for 
the interacting theory in Fig. || for jx = 0.5. The left figure shows the Langevin history of (|0| 2 ), 
using three different initial conditions but the same random number sequence in the stochastic 
process. Since /I = 0.5 is far from the critical value (and hence the minimal value of A p is well 
separated from zero), the memory of initial conditions is quickly erased. To have an indication for 
the size of the statistical fluctuations during the evolution, runs with a different random number 
sequence are shown in the figure on the right. Note that in a typical run the Langevin time < 6 < 
250. 




1 2 3 4 1 2 3 



Langevin time Langevin time 



Figure 2: Independence of initial conditions: Langevin evolution of (|0 1 2 ) with three different initial condi- 
tions using the same (left) and a different (right) random number sequence (ji — 0.5, m = A = 1, N = 4 4 ). 

Expectation values are obtained by noise averaging and taking the limit 8 — > °°. Provided 
again that A p > 0, we find from the analytical solution 

1 2A 2 -\- B 2 

Km (^(0)^(0)) = S ab S pi/ — 1 ^f, 

iimjCp(e)^(e)) = ^V^ I i^i, 

hmj$_ p ml p ,(d)) = e ab 8 p /--^^. (2.7) 
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Figure 3: Comparison between the mean field predictions (lines) and the numerical data (symbols) for 
| 2 ) in the full (left) and phase quenched (right) theories in the Silver Blaze region. The vertical dotted 
line indicates the mean field estimate for the critical chemical potential. 



So far interactions have been neglected in this analytical calculation. An elegant way to include 
interactions on the mean field level is by looking for fixed points of the Langevin equations for the 



two-point functions above, using a Gaussian factorization [10]. After the dust settles, the result is 
a simple shift in the effective mass parameter, m 2 — » m 2 + 4A(|0| 2 ), as expected. The critical value 
of the chemical potential in the mean field approximation then follows from m 2 + 4 A ( | <j) \ 2 ) + 2 ( 1 — 
coshjit) = 0, where (\<p\ 2 ) is determined by a self-consistent gap equation. For the parameters used 
here we find }i. c = 1.15, in agreement with the numerical estimate. 

We can now compare the nonperturbative results from the numerical solution of the complex 
Langevin equations with the mean field estimates on a finite lattice. This is shown in Fig. ^] for 
(|0| 2 ) and in Fig. |] for (n). Since at /I = the action is real, the value of (|0| 2 ) at /I = is 
obtained with real Langevin dynamics. This provides a nontrivial check of the numerical code. We 
observe excellent agreement between the simulations and the analytical results, including finite- 
size effects. This indicates that complex Langevin dynamics is successful, but also that the mean 




Figure 4: As above, for the density (n) 
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Figure 5: Left: difference Af between the free energy densities of the full and the phase quenched theories 
in the mean field approximation. Right: as in the preceding figures, for the average phase factor (e'^pq in 
the phase quenched theory. 

field approximation is adequate and the theory is effectively weakly coupled. 

The results so far answer three of the questions raised in the Introduction. What about the 
severity of the sign problem? From the pronounced difference between the full and the phase 
quenched theory, one can expect the sign problem to be severe. This is quantified by studying the 
average phase factor, 

ei<? = ]^y (2 " 8) 

where the expectation value is taken in the phase quenched theory, with real weight \e~ s \. This 
phase factor is relevant for reweighting and is expected to vanish exponentially in the thermody- 
namic limit, since 

<^ )pq = |* =e -<W (2.9) 
z pq 

where Af is the difference between the free energy densities in the full and the phase quenched 
theories. The mean field results for the difference between the free energies are shown in Fig. || 
(left). We note a rapid convergence in the thermodynamic limit. Inserting these results for Af in the 
exponential e~ a ^ yields the average phase factor. In the numerical simulations the phase factor 
is measured directly. A comparison between the two is shown in Fig. || (right). We observe again 
excellent agreement. The severeness of the sign problem is manifest. As mentioned above, this is 
necessary and does not affect complex Langevin dynamics in practice. 

3. Outlook 

The question posed in the title and made more specific in the Introduction has been answered 
positively for the relativistic Bose gas at finite chemical potential considered here. In this model 
the sign problem is required to yield the correct physics, summarized by the Silver Blaze problem, 
and it is severe. Nevertheless, complex Langevin dynamics can be applied without encountering 
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(numerical) problems affecting e.g. stability or convergence. The outcome of the numerical simu- 
lations can be understood from a comparison with mean field theory. Moreover, the inner workings 
of the algorithm can be analysed as well, relying essentially on the weakly coupled nature of this 
theory. 

Does this mean that stochastic quantization can easily be applied to other theories with a 
complex action due to a nonvanishing chemical potential? We are currently studying complex 
Langevin dynamics in a variety of models, including the three-dimensional XY model at nonzero 
chemical potential, and extending our analysis of QCD with static quarks. At this stage it is fair to 
say that easily is slightly too optimistic. 
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